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Abstract. Kinetic equilibration during the early stage of a relativistic heavy 
ion collision is studied using a radiative transport model. Thermalization is 
found to dominate over expansion with medium regulated cross sections. Pres- 
sure anisotropy shows an approximate as scaling when radiative processes are 
included. It approaches an asymptotic time evolution on a time scale of 1 to 2 
fm/c. Energy density is also found to approach an asymptotic time evolution 
^ ' that decreases slower than the ideal hydro evolution. These observations indi- 

^ , cate that viscosity is important during the early longitudinal expansion phase 

of a relativistic heavy ion collision. 
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O ' 1. Introduction 

0^ ■ 

' One of the evidences for the discovery of the Quark-Gluon Plasma is the collective 

behavior of produced particles that is consistent with ideal hydro description of the 
evolution of the hot and dense matter produced in relativistic heavy ion collisions. 
^ . , To go beyond the ideal hydrodynamic description, viscous hydro has been developed 

■ [ [11 [21 [3l m [5l H] . From a different perspective, microscopic transport models 

have also been used to study the degree of thermahzation [[71 [3 [9]. These studies 
show that even small viscosity can affect the evolution and can be important for the 
precise interpretation of experimental data. In the following, we will describe kinetic 
equilibration of a Quark-Gluon Plasma. In particular, we will discuss what happens 
when medium regulated cross sections are used and when chemical equilibration via 
particle number changing processes is included. We will look at the central cell in 
central heavy ion collisions. The kinetic equilibration will be characterized by the 
pressure anisotropy. In addition, the energy density evolution will also be studied. 
Finally, a summary will be given. 
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2. Pressure anisotropy 

Local thermal equilibrium requires local pressure isotropy. For the central cell in 
central heavy ion collisions, the degree of pressure isotropization can be character- 
ized by the longitudinal pressure to transverse pressure ratio. Here the longitudinal 
direction is the beam direction. Because of azimuthal symmetry, the two transverse 
directions have the same pressure and the average will be used as the transverse 
pressure. When the system is isotropic, the longitudinal pressure is expected to be 
the same as the transverse pressure, and the pressure anisotropy is expected to be 
equal to one. We will first review findings from the AMPT model, then we will 
introduce the new ingredients for the thermalization study, i.e., medium regulated 
cross sections, and particle number changing processes. Then we will look at how 
the two new ingredients can help to maintain kinetic equilibrium by starting with 
local thermal initial conditions. Furthermore, we will start with zero initial pressure 
anisotropy as produced from an inside-outside parton production process. This will 
give information on how much equilibration a model can generate with the inclusion 
of these two new ingredients. 

Kinetic equilibration from the AMPT model was studied in Ref. [fSl. The 
AMPT model is a hybrid model that uses different models for different stages of 
relativistic heavy ion collisions [[TUl [HI HH [13]. It has two versions: the default 
model and the string melting model. The default model evolves minijets in HIJING 
[ [T4] with ZPC [[15] and reconnects them with their parent strings for fragmentation 
with PYTHIA [[Ej. The string melting model dissociates hadrons from HIJING 
into their constituent quarks and anti-quarks. After the evolution of the quark-anti- 
quark plasma with ZPC, a coalescence model is used to turn partons into hadrons. 
The hadron stage is handled with the ART [ [TT] hadron cascade. The AMPT 
model has been used to study a variety of experimental observables [[18] [TH [201I2T] . 
The kinetic equilibration study yields a few interesting features. First, the string 
melting model has a more rapid increase of the pressure anisotropy with proper 
time as compared to the default model. The larger the cross section is, the faster 
the increase of the pressure anisotropy. The pressure anisotropy reaches a plateau. 
Later, as transverse expansion sets in, it increases and crosses one. It does not stay 
at one for any significant time period. This shows that only partial thermalization 
is achieved in these collisions. 

The AMPT model has fixed cross sections for partons. In the following, we will 
make use of medium dependent cross sections and add particle number changing 
processes. The two gluons to two gluons (2 to 2) cross section is taken to be the 
pQCD cross section regulated by a medium generated Debye screening mass. It is 
related to the screening mass by 



1722 = 



(1) 



The 2 to 2 cross section (722 is proportional to the strong coupling constant as 
squared and inversely proportional to the screening mass fi squared. The screening 



Kinetic equilibration 



3 



mass is related to as, the cell volume V, and the particle momenta via 

Notice that as the number of particles increases, even with the same momentum 
distribution, fj? increases. This leads to a smaller 2 to 2 cross section. In other 
words, the collision rate is regulated by screening, particles tend to interact with 
close particles and can not reach very faraway particles. This is a desirable feature 
and can avoid conceptual and technical problems related to large cross sections in 
dense medium. 

The simplest particle number changing processes will be included. These are 
the 2 to 3 and 3 to 2 processes. The 2 to 3 cross section is taken to be 50% of the 
2 to 2 cross section. This is roughly the ratio obtained in a more elaborate model 
with pQCD matrix elements [[22]. The 3 to 2 rate is determined by detailed balance 
to ensure correct chemical equilibration. The 3 to 2 reaction integral is related to 
the 2 to 3 cross section by 

/32 = 192^2^23. (3) 

This relation assumes isotropic reactions. In the following, the outgoing particle 
distribution will be taken as isotropic. This is expected to catch the gross features 
of particle interactions in the very dense central region while for the peripheral 
region this may lead to deviations from results obtained from more realistic forward 
type particle distributions. 

We will start with local thermal initial conditions at an initial proper time 
T — 0.5 fm/c. In particular, two initial temperatures (Tq) will be used, 0.5 GeV and 

1 GeV. The space-time rapidity will start uniformly within a range from -5 to +5. 
The initial transverse coordinates will be within a disk of radius of 5 fm. The space- 
time rapidity density is taken to be dN/drjs = 1000. The low initial temperature 
case is similar to what is expected at RHIC while the high temperature case is a 
sketch of what may happen at the LHC. (Note that the space-time rapidity density 
increases toward 2000 in the high temperature case with its large initial energy 
density.) For the early evolution, a time period from to 2 fm/c will be used. This 
is early enough and no transverse expansion is expected to affect the results. The 
evolution in the central cell can be considered as having only longitudinal expansion. 

Fig. [T] shows the proper time evolution of the pressure anisotropy with the 
above set up. The free streaming curve can be used as a guide of the effect of 
longitudinal expansion. It goes slightly slower than l/r^. As the interactions are 
turned on, the decrease of pressure anisotropy slows down and after some time, the 
pressure anisotropy reaches a minimum and then it increases with time trying to 
return to isotropy again. As the coupling constant increases from 0.3 (the weak 
coupling case) to 0.6 (the strong coupling case), the pressure anisotropy goes closer 
to 1 and kinetic equilibration becomes faster which can be seen from the fact that 
the minimum happens earlier. A closer look at the evolution shows that if only the 

2 to 2 process is included, there is an a^To scaling [ 23 . If one evolution is solved, 
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Fig. 1. Pressure anisotropy evolution ^ig. 2. Like Fig. [T] with additional cases 
from local thermal initial conditions with t^j^t have a higher initial temperature, 
different coupling constants. 



then the evolutions of systems with the same agTo are all solved. As 2 to 3 and 
3 to 2 processes are turned on, the evolution instead follows an approximate as 
scaling (Fig. [2|). For the same as, the higher the initial temperature, the closer the 
pressure anisotropy is to one. However, as a result of the smaller cross sections due 
to particle production, the high initial temperature weak coupling curve does not 
go high enough to match the low initial temperature strong coupling case as in the 
case with only elastic collisions. 

Next we will look at what happens when there is no fast thermalization. We 
will change the initial particle anisotropy. At the initial time, the directions of 
particle momenta will be changed so that they all move in the transverse plane. 
This is consistent with the inside-outside particle production picture. It can be 
compared with the case with an isotropic initial particle distribution and the same 
initial energy density. Fig. [3] shows that the two cases with two different initial 
pressure anisotropics but the same initial energy density move toward a common 
asymptotic evolution. In other words, as the memory of the initial condition fades 
away, the system evolves toward some common evolution. The memory effect can 
be measured by the difference between these two curves at any proper time relative 
to the initial difference. A relaxation time can be used to describe how fast the 
memory gets lost or how fast the isotropization is. We notice that as the coupling 
constant increases, the relaxation time decreases. Likewise, as the initial energy 
density increases, the relaxation time decreases. Next we will see if some common 
evolution also happens in the energy density evolution. 
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X (fm/c) (f^i/c) 

Fig. 3. Like Fig. El with additional cases ^ig- 4. Energy density evolution on the 
that start with zero pressure anisotropy. log'log scale. 



3. Energy density evolution 

For the evolution of azimuthally symmetric system of massless particles, there are 
two independent thermodynamic quantities for the central cell. In addition to the 
pressure anisotropy, the energy density can be chosen as the other variable. 

Fig. [5] is a log-log plot of the energy density evolution. We will get a closer 
look at the case starting with the high initial energy density. The low initial energy 
density case shows roughly the same behavior. First, we notice that the curves 
starting from the same initial energy density are bounded by two curves, the upper 
dotted one is the free streaming evolution from the transverse initial condition. The 
lower solid one is the isotropic initial condition with strong coupling. If particles free 
stream from the transverse initial condition, they will stay within their cell when the 
system expanses longitudinally. The total energy in the cell remains the same while 
the volume increases due to longitudinal expansion. The volume is proportional to 
the proper time. So the energy density decreases as l/r. The evolution from an 
isotropic initial condition with strong coupling constant can be fitted by a power 
law. It gives e ex l/r^-^^. 

The lower dotted curve is for free streaming from the isotropic initial condition. 
It follows the isotropic initial plus strong coupling curve in the beginning and devi- 
ates from it to approach the same power law behavior for free streaming from the 
transverse initial condition. Other curves also follow the early time power law be- 
havior which reflects the initial anisotropy. They also develop asymptotic power law 
behavior that reflects interactions at late times. The upper solid and upper dashed 
curves are from the transverse initial condition with strong and weak interactions 
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respectively. They fohow 1/t at early thnes and approach l/r^-^^ at late times. 
The strong coupling case has smaller energy density and the evolution is closer to 
the strong coupling case from the isotropic initial condition. The lower solid and 
dashed curves are for evolutions from the isotropic initial condition with strong and 
weak interactions respectively. Both follow l/r^-^^ approximately while the strong 
coupling case has slightly smaller energy density at late times. The time period 
that we studied is about the time period when there is only longitudinal expansion 
for our system. As from the pressure anisotropy study, this is not much longer than 
the relaxation time. As a consequence, the time evolution is slightly slower 

than the ideal hydro l/r^ '^^ energy density evolution. Hence, viscous corrections 
are present and important during this time period. 

4. Summary 

In this study, pressure anisotropy is shown to result from the competition between 
longitudinal expansion and kinetic equilibration. With medium regulated cross sec- 
tions, thermalization dominates over expansion at late times. There is an approx- 
imate as scaling for the evolution of pressure anisotropy when inelastic processes 
are included. When evolving from the same initial energy density, different initial 
pressure anisotropics follow the same asymptotic evolution when memory effect goes 
away. The relaxation time is on the order of 1 to 2 fm/c. At the same time, energy 
density also approaches the limiting power law evolution. This evolution is slower 
than the ideal hydro evolution. This indicates the importance of viscosity during 
the early longitudinal expansion stage of evolution. Further studies with more so- 
phisticated initial conditions and interactions will help to improve the estimate of 
viscous effects and the impacts on experimental observables. 
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